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Abstract 



A phenomenological analysis of the general two Higgs doublet model is 
presented. Possible constraints of the Yukawa couplings are obtained from 
(-^ : the K^ - K^, B^ - B^ and D^ - D° mixings. A much larger D^ - D° 

!>■ , mixing than the standard model prediction is possible. It is shown that the 

emerging of various new sources of CP violation in the model could strongly 
cn '_ affect the determination of the unitarity triangle. It can be helpful to look 

for a signal of new physics by comparing the extracted angle /? from two 



Q^ ■ different ways, such as from the process B — > J/ipKs and from fitting the 

5j ■ quantities \Vub\, Aitib and e. 

P^; PACS numbers: ll.30.Er, 12.60.Fr 
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I. INTRODUCTION 

In the standard model(SM) of an electroweak SU{2)i x U{1)y gauge theory with 
only one Higgs doublet, the only source of CP violation comes from the complex Yukawa 
coupling between Higgs and fermion fields ^. Since the Higgs sector of SM is not well 
understood yet, many possible extensions of SM have been proposed 0. One of the 
simplest extensions of the SM is to simply add one Higgs doublet without imposing the 
ad hoc discrete symmetries. For the convenience of mention in our following discussions, 
we may call such a minimal extension of the standard model with generally adding an 
extra Higgs doublet as an S2HDM p-pT| and assume CP violation solely originating from 
the Higgs potential [T^J7|,3]- The most general Yukawa coupling and Higgs potential can 
be written as: 

Ly = gL(rf 01 + r2^02)t/ij + gL(rf 01 + t^<P2)Dr (i) 

and: 

1^(01, 02) = - ^?0101 - ^A.4>2<P2 - (^?20l02 + h.c) 

+ Ai(0l0i)2 + A2(0502)' + A3(0l010502) + A4(0l020^20l) 



+ ^[A5(0l02)' + h.c] + [(A60l0i + A70502)(0l02) + h.c] (2) 

The major issue with respect to the two Higgs doublet model is that it allows flavor 
changing neutral current (FCNC) at tree level, which must be strongly suppressed in 
K^ — K^ and B^ — B^ mixing processes. In order to prevent FCNC from tree level, an 
ad hoc discrete symmetry is often imposed 

01 —* -01 and 02 -^ 02 
Ur, -^ -Ur^ and Dr^ -^ tDr^ (3) 

Thus, one obtains the so called mode I and model II, which depends on whether the up- 
type and down-type quarks are coupled to the same or different Higgs doublet respectively 
0. Once the discrete symmetry is adopted, the factor fiu, Xe and A7 in Eq.(2) mush 
vanish. As a result, no CP violation can occur from V{(f)). Thus the only source of 
CP violation is the complex Yukawa couplings, which lead to a phase in the Cabibbo- 
Kobayashi- Maskawa(CKM) quark mixing matrix. 

In contrast, one can replace the discrete symmetry by an approximate global family 
symmetry [§,||,|^J^ , thus the suppression of FCNC can be explained via the smallness 
of the off-diagonal terms. Furthermore, when abondoning the discrete symmetries, one 
can obtain, after spontaneous symmetry breaking, rich sources of CP violation from a 
single relative phase between the two vacuum expectation values of Higgs fields. It has 
been shown |j^J^ that even when the CKM matrix is real, the single phase arising from 
the spontaneous symmetry breaking can provide enough CP violation to meet the exper- 
imental measurments. One of particular important observations is a new source of CP 
violation in charged Higgs boson interactions, which is independent of the CKM phase 
and can lead to a value of e'/e as large as 10^'^ [Q. In the S2HDM, the two Higgs fields 
have, in general, the vacuum expectation values: 



< 05 > = ^ COS /?e*' 
v2 



aO ^ 



<0^> =^sin/3 (4) 



It is natural to use a suitable basis: 



Hi = cos/?0ie *^ + sin/302 

H2 = sin PM''^ - cos /502 (5) 

such that: 

where H'^,R,I are real Higgs bosons. The three neutral scalars H^ = {R,p,I) can be 
rotated into mass engenstates hi = {h, if °, A) via orthogonal matrix O^ . 

Hi = OJf,Hl (7) 

From approximate global family symmetries, we know that the Yukawa coupling ma- 
trices rf in Eq.(l) have small off-diagonal elements, typically between 0.01 and 0.2 in 
order to meet the constraint of FCNC from K^ — K^,B^ — B^ mixing. The Yukawa 
interaction can be rewritten as [0: 



LY = {Li + L2){V2GFf 



/2 



with 



Li = V2{H+J2^^'i,Vi,uld^-H-J2^^^m^^Vidiui 






3 

i 

3 3 

+ {R + il) Y. UmdAd^R +{R- il) E ^u^'^u^u'r + h.c (9) 

i i 

(3 3 

+ {R + il) E 4-44 + {R- il) E ^^J^L^R + h.c, (10) 

where Li has no flavor-changing effects other than that expected for H^ from the CKM 
matrix V and L2 contains the flavor- changing effects for neutral bosons as well as small 
additional flavor- changing terms for if^. The factor ^f^mf. and /if- arise primarily from 

the diagonal and off-diagonal elements of T{ respectively. 



There are four major sources of CP violation [0,|1: (1) CKM matrix,(2) the phases in 
factors ^f. which provide CP violation in charged- Higgs boson exchange; (3) the phases 
in fx^j, this yields CP violation in FCNC. (4) CP violation in mixing matrix O^. One 
of the most distinctive features of these sources is that the factors ^j. can provide CP 
violation in charged Higgs boson exchange in addition to and independent of CKM phase. 
As a consequence, in AS* = 1 transitions its contribution to e'/e could be as large as 
10^'^ and become comparable with the experimental data [0,in|. Thus a measurement of 



e'/e would not necessarily be due to CKM mechanism |7|J^,p!5|l. This paper is organized as 
follows: In section 2, the analysis of the possible constraints of the Yukawa couplings from 
F^ — F° mixings is presented. The various sources of CP violation and their influence on 
the determination of unitarity triangle are discussed in section 3, and section 4 contains 
our conclusions. 



II. CONSTRAINTS FROM K^ - K^, B^ - B^ AND D^ - D^ MIXINGS 

In the standard model, it is known that the neutral meson mixings arise from the box 
diagram through two-W-boson exchange. The extremely small values of the neutral K and 
B mass differences impose severe constraints on new physics beyond the SM, especially 
on those with FCNC at tree level. In the S2HDM, additional contributions to the neutral 
meson mixings can arise from the box diagrams with charged-scalar exchanges and tree 
diagrams with neutral-scalar exchanges. The mass difference of K^ — Ks is given by 

AniK ^ 2i?eMi2 = 2Re{M^^ + Mg" + Mg^ + Mg° + M[^) (11) 

where M]^^, Mg^ and M^^ are the contributions from box diagrams through two 
W^- boson, two charged-scalar H^ and one W- boson and one charged-scalar exchanges 
respectively. Mf^ is the one from the FCNC through neutral scalar exchanges at tree level 
. M[2 presents other possible contributions, such as two-coupled penguin diagrams and 
nonperturbative effects. They are resulted from the corresponding effective Hamiltonian 

i^eTf = -^2'^w E %A,A,V5^i?'^^(^., x,)dj,{l - i,)sdY{l - 75)5 (12) 

■f^7/.(l - 75)s^7''(l - 75)5 + ^/XsXdy/y7ydilifdlt,il + 75)s^7''(l + 75)5 
+25ijy/x~x'j^/y7y~disCdi^Cjdl^l{.^ + 75)5^7^(1 - 75)5] (13) 

+5^^(y„y,)Vx^[a;,ef e;e;rfli - i^)sd{i - ^,)s 

+Xs^%^jdil + 75)sd(l + 75)5 + 2./^E:rdU*d^iCjdil + 75)srf(l - 75)5]} 

id 

■dj.ii - 75)5^7^(1 - 75)5 + iy^ + yJ)V^:I^sU*d[B?'^{y^, y„ yu,) (14) 

■da^,{l - i,)sda^''{l + 75)5 + 5^^(y„ y„ y^)d{l - ^,)sd{l + 75)5]} 



where the B^^ , Ef ^, B§^, B^^ , B§^ and B^^ arise from the loop integrals, and 
they are the functions of Xi = mf/m^r and yi = mf/m\j with i = u,c, t, W, their explicit 
expressions are presented in the Appendix. Here r]ij, r]fj^ and r]fj^ are the possible QCD 
corrections and Aj = VisV*i. Note that in obtaining above results the external momentum 
of the d- and s-quark has been neglected. Except this approximation which should be 
reliable as their current masses are small, we need to keep all the terms. This is because 
all the couplings Aj and C,i are complex in the model, even if some terms are small, they 
can still play an important role on CP violation since the observed CP-violating effect in 
kaon decays is of order 10~^. The contributions of neutral Higgs bosons exchange at tree 
level can be evaluated by 

M»" = < i'»ii/f;;ip» > 



with 



127r2 y mj, m^. ■'j ^ mf^omf^ 



1 i 



Sk,ij is related to /if,- through: 



sL = {0^k + ^^f0^k)^= (16) 



/ 



mjm 



J 



where o"/ = 1 for d-type quarks, o"/ = — 1 for u-type quarks. The formula is expressed 
in a form which is convenient in comparison with the one obtained from the box diagram 
in the SM, here Jrrifjmf. with convention i < j plays the role of the CKM matrix element 
Vij, and mj' is introduced to correspond to the loop-quark mass of box diagram. Namely 
/' and fj are the two quarks in the same weak isospin doublet. Note that the result is 
actually independent of nifi . Here nif. are understood as the current quark masses. In 
our following numerical estimations we will use niu = 5.5 MeV, m^^ = 9 MeV, m^ = 180 
MeV, rric = 1.4 GeV and rrib = 6 GeV with being defined at a renormalization scale of 
1 GeV. fpo and mpo are the leptonic decay constant (with normalization /^ = 133MeV) 
and the mass of the meson P^ respectively. Bpo and fpo are bag parameters defined via 

< p'mi±i.)mp' >= - . ^^':^^\, ^po (17) 



<P0|/.(l±75)/,/.(l±75)/,|P°> ^ ^ 

In the vacuum saturation and factorization approximation with the limit of a large number 
of colors, we have Bpo -^ 1 and fpo -^ 0, thus Y/j^^j = Zl^^ 

It is known that H^f^ contribution to ArriK is dominated by the c-quark exchange and 
its value is still uncertain due to the large uncertainties of the hadronic matrix element 



< K''\{d^^{l-^,)sy\K' >= --f^mj^BK (19) 



where Bk ranges from 1/3 |]T6[ (by the PCAC and SU{3) symmetry), 3/4 [T^ (in the 



hmit of a large number of colors) to 1 [jT9[ (by the vacuum insertion approximation). The 
results from QCD sum rule and Lattice calculations lie in this range. For small Bx, the 
short- distance H^^ contribution to Auik fails badly to account for the measured mass 
difference. 

In general, we obtain 

AfUK = p^fxBKmKmlsm'9{r^,,B'^'^{x,) + ^T^^/y^l^^l'B^^'iy,) 
bT[ 4 



which is subject to the experimental constraint [|1? 



ry2 

ArriK = 3.5 x IQ-^el/ ~ v^—^/im^m^ sin^ 6 (21) 

The effective Hamiltonian for B'^—B^ Mixing is calculated with the aid of the box diagrams 
in full analogy to the treatment of the K'^ — K^ system. Its explicit expression can be 
simply read off from the one for K^ — K^ system by a corresponding replacement: s ^->- b. 
The "standard approximation" made there, namely neglecting the external momenta of 
the quarks, is also reliable since dominant contributions come from the intermediate top 
quark. With this analogy, the considerations and discussions on K'^ — A'° mixing can 
be applied to the B^ — B^ mixing for the contributions from box diagrams. As it is 
expected that |ri2|/2 <C |Mi2| in the B-system (which is different from K-system), the 
mass difference for B'^ — B^ system is given by Ams ~ 2 1 M12 1 . 

The general form for the mass difference in the B^ — B^ system can be written 

r'2 , 1 1 



IS ^ ^^{fB^BBriuYmBmi\Vu?-\{riuB'^'^{xt) + -nf,^ y^W B^"" {y,) (22) 

Gtt^ " rjtt 4 

+2% yt\it\By (yf,yj] + — )^( ) —772(^^,13) } I 



which is subject to the experimental constraint [^ 



ry2 

AniB = (3.1 ± 0.12) X lO^^eF ^{125MeVy niBinQGeVf {sin 9 = 0.22)^ (23) 

It is known that in the standard model the short-distance contribution to ArriD from 
the box diagram with W-boson exchange is of order of magnitude Am^°^ ~ 0(10^^) eV, 
here the external momentum effects have to be considered and were found to suppress 
the contribution by two orders of magnitude ^^. This is because of the low mass of 
the intermediate state. It is not difficult to see that the additional box diagram with 
charged- scalar exchange gives even smaller contribution except \^s\ is as large as |.^s| ~ 
2mH+/rns which is unreliable large for the present bound mH+ > 54.5 GeV ||^. It 
has been shown that dominant contribution to Am£, may come from the long-distance 
effect since the intermediate states in the box diagram are d- and s-quraks. The original 



estimations were found that Arrio ~ 3 x 10 ^ eV [^ and ArriD ~ 1 x 10 ^ eV 



An alternative calculation ||2J] using the heavy quark effective theory showed that large 
cancellations among the intermediate states may occur so that the long-distance standard 
model contribution to Arrio is only larger by about one order of magnitude than the short- 
distance contribution, which was also supported in a subsequent calculation [p5| . 

With this in mind, we now consider the contribution to Atud from the neutral scalar 
interaction in the S2HDM. It is easy to read off from eq.(15) 

Am, = 2|M,| = ^JnBnmni^-) m,Y.{———) \Y,,,\ 

= 0.64 X io-^( W|Z)2 '500GeV,\Y^ 

With the above expected values in the second line for various parameters, the predicted 
value for Amo can be closed to the current experimental limit |Amj:)| < 1.58 x 10~^ eV 
|T^, this implies that a big D° — -D° mixing which is larger than the standard model 



prediction is not excluded. With this analysis, we come to the conclusion that a positive 
signal of neutral D meson mixing from the future experiments at Fermilab, CESR at 
Cornell and at a r-charm factory would be in favor of the S2HDM especially when the 
exotic neutral scalars are not so heavy. 

We now proceed to study the constraints on the parameters of the model. Since the 
parameters ^j- and fj.{j are in general all free parameters, for simplicity we will consider 
the constraints in two extreme cases. 

Case 1: the mass difference is purely explained through the additional box diagrams 
from two scalar-boson and one W-boson one scalar-boson. In this case, the parameter 
^/. is of particular importance. Both its amplitude and phase will play an important 
role in the neutral meson mass difference and CP violation. It is quite different from 
the earlier analysis in type 1 and type 2 2HDM p6[ in which the three couplings C,u, 6 



CI 



and ^t are equal, i.e., ^u = ^c = ^t = tan/3 . This is why the constraint from e is much 
stronger than the one from Arrix in those models. In the S2HDM, one has in general 
iu 1^ ic ^ iti there are more degrees of freedom to fit e and Am^, ArriB as well as Amu 
. Since the main contribution to Attlk comes from c-quark though the loop, the upper 
bound of ^c can be extracted from K^ — K^ mixing. The result is plotted in Fig.l. In the 
numerical calculations we take fx = 161 MeV and Bk = 0.75. The rang of m^ is from 
100 to 500 GeV. Since the bound of |^c| strongly depends on the SM prediction on Amx, 
three different values oi Aqcd{-^qcd = 0.21, 0.31, 0.41) have been used, the corresponding 



ratios to the experimental data {AmK)exp are 0.52,0.67 and 0.91 |2^ 

In the B system. It is of interest to study its relative ratio to the SM, since large degree 

of uncertainty can be avoided from CKM matrix IV^^^I and hardonic matrix elements. In 

Fig. 2 we illustrate the relation between \C,t\ and the charged Higgs mass m'jj when the 

ratios of the HW and HH box diagram contributions to the WW-box diagram contribution 

in the SM are 2:1, 1:1 and 0.5:1. 

As it is shown in Fig.l and Fig. 2, in general \C,c\ is much larger than |,^i|. Even when 

the HW and HH contributions to ArriB are twice as much as the SM one, \^c\ can still be 

larger than |^i| by an order of magnitude. 
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FIGURES 
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FIG. 1. The upper bound of |^c| with respect to the mass of the charged Higgs scalar in 
case 1. The three curves are corresponding to the ratios {/^mK)sM /{^fnK)exp = 0.52 (dotted), 
0.67 (dashed) 0.91 (sohd) 



M„HGeV} 
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FIG. 2. The value of |.^t| with respect to the mass of charged Higgs m^. The three curves 
are corresponding to the different ratios of HW and HH box diagram to the one from the SM: 
2:1 (dotted), 1:1 (dashed), 0.5:1 (solid). 



Case 2: the mass difference is fitted through neutral-scalar exchange at tree level. 

If S-j is 
expected to be symmetric under the exchange i '-^ j and rpo = 1, F/j • gets the following 



From Eq.(15) we know that the parameters arise in M^2 ^^^ ^kij rather then /ifj. 



simple form 



yf 

k,t,j 









+ ^fO!^k 



_Ref4j_ 



(25) 



Hence both imaginary and real part of /if^ are of importance. Futher more, the phases 
in nlj also provide a new source of CP violation as we have mentioned in the previous 
section. To simplify the discussions, we assume that one of the scalar-bosons, for example, 
the scalar h is much lighter than the other two scalars H and A. The scalar bosons H 
and A are assumed to be heavier than 500 GeV. The upper bounds can be obtained from 
j{0 _ j^o^ j^o _ j^o g^j^j^ jjo _ ^0 mixings. The present consideration is more general than 
the one in |P] where all the couplings Yl^ij are setted to be equal. As a consequence, the 
constraints from different meson mixings provide different upper bounds upon different 
y/j-. The results are shown in Fig.3. It is seen from the fig.3 that the upper bound of 
Y^-^2 is much higher than that from K^ and B^ system. This implies that a larger D^ — D^ 
mixing than the standard model prediction is possible. 



100 150 200 250 300 350 400 450 
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M^t(GeV) 



FIG. 3. The m^ dependence of the upper bound of ReYj^--. ReYfi2 from K^ — ^^(solid), 
ReVfi^ from B^ — ^'^ (dashed) and ReY^i2 from D^ — Z)*^ (dotted). The mass of the other scalar 
mP^ is fixed at m\ = 500 GeV. 



III. CP VIOLATION AND UNITARITY TRIANGLE 



Besides the neutral meson mass difference, the indirect CP violation parameter ex 
could also provide constraints on the values of ^/. and /if^. The standard definition of e is 



ImMi2 . , ^ j^/A 



v^'2i?eM: 



+ io)e'-" (26) 



12 



where ^o = ImAo/ReAo with \Ao\ = (3.314±0.004) x 10^''' GeV the isospin-zero amphtude 
oi K ^ TTTT decay. Usually, the ^o term is relatively small as it is proportional to the small 
direct CP- violating parameter e'. 

The first part of contribution to e comes from the box diagram through W-boson and 
charged-scalar exchanges 

f^2 c,t 



f'^BKiriKmirrijlY^ Im{XiXj)ReBij{mi, mf, ^i, ^j 



12vr2 

+Re{\i\j)ImBij{mi, rrij; ^i, Q} (27) 

where Bij{mi,'mj]^i,^j) depend on the integral functions of the box diagrams and their 
general form is given in the Appendix. The imaginary part ImBij{mi,mj]C,i,C,j) arises 
from the complex couplings ^j. 

The second part is due to the flavor changing neutral scalar inteactions at tree level 

/™M- ^ ^/iB,,,„,.(g)^,„J S(?|^)^/™(>'i.)' (28) 

This provides a contribution to e in almost any models which possess CP-violating flavor 
changing neutral scalar inteactions. 

In particular, the parameter e could receive large contributions from the long-distance 
dispersive effects through the vr, r] and rj' poles [28|. For a quantitative estimate of these 
effects, we follow the analyses in refs. p8|-pl[l 



TT,r],r) 



Im{< K^\L^ff\i >< i\Leff\K^ >) 



1 2k 



<ir°|L_|7r°><7r°|L+|7r''>) (29) 



ArriK mj^ — m^ 

G^ 2 D/ /"^i^N2 • a 2/ I^Ois 3KAKn 



,f],B'^m^C^r Sin 6ml 



■ ^[/mAii?ei^^(mi, ^0 + ReXJmP,''{mi, ^0] (30) 

i 

where k is found to be k ^ 0.15 when considering the SU{3)— breaking effects in the 
K—r]s transition, and nonet-symmtry- breaking in K—rjo as well as 77—77' mixing. We shall 
not repeat these analyses, and the reader who is interested in it is refered to the paper 
PT| and references therein. L_ and L+ are CP-odd and CP-even lagrangians respectively 
(with convention Lg// = L+ + iLJ). The L_ is induced from the gluon-penguin diagram 
with charged-scalar 

L- = fsda^,{l + l,)X''sGl, - hda^,{l - 75)A"sG;, (31) 

with 

10 



f, = ^^2^^ E Im{U,X.)y,P?{y,) (32) 

where P^iiji) is the integral function and presented in the Appendix. From fs and fd 
it is not difficult to read off the ReP^^ {mi,^i) and ImP^irrii^^i) (see Appendix). In 
obtaining the last expression of the above equation, we have used the result for the 
hadronic matrix element < K^\L_\'k^ >= (/^ — fdl^K-K where Ak-k has been computed 
in the MIT bag model and was found |3^ to be Ak-k = 0.46*61^^ for a^ = 1, and the 
convention < 7r°|L_|_|_ft'° >= |G/|^i?^m|-(2m7^/ms)^ sin6', where 5^ is introduced to fit 
the experimental value < 7r°|L_,_|K° >= 2.58 x 10^'^GeV'^ and is found to be S^ = 1.08. 
We then obtain y/Wa^SnAKw / [4^V^'m'si'm'K — w^^)] — 1-4. 

Neglecting the t— and m— quark contributions and also the terms proportional to rrid in 
comparison with the terms proportional to rris, the total contributions to the CP- violating 
parameter e can be simply calaculated from the following formula 

|6| = 3.2 X 10-^5;.(^)^^||^sin5^M{-^[r/cc5'^^(xe) 

+\v^/yc\U'B^''iyc) + 2vryc\U'Briyc,y^)] 

, AVcb\fnt 2f-. \Vub\ X NT Tyww, 



2mc \Vcb\\Vus\ 

+\vrytUB^''iyt) + 2vryt\^^B^'^{yuy^)] 
+^^[vctB^^{xc,x,) + lv^,''V^t\^c\%\'B^''{yc,yt) 

+^vTVy^tRe{Ut)B^'^{yc, yu y^)]} (33) 

+2.27 X 10"^ ^ ^'^^' BKYi Y 

6.4 X 10-3 ^ V mHO ' 

^ k 

+2.27 X lQ-'lm{CCf-^^BK{ln'^ - 2) 

+2.27 X 10-3/m(e.e.)^^i?K(^n^ - |) + 4 

mjj+ mf. 2 V2 

where we have used the experimental constraint on 2ReMi2 = Am^^' 

It is analogous to the section 2, we consider the contributions to e in two different 
cases. The ffist one, CP violation is governed by the induced KM mechanism, i.e., the 
ffist term of the above equation becomes dominant. In this case, new contributions come 
from the box diagrams of two charged-scalar, and one W-boson and one charged scalar 
exchanges. Since the expression contains Re{^c^t), the relative phase 9 between ^c and 
^t, i.e., Re{C,cC,t) = \(,c\\C,t\ cos9, may play an important role. It is of interest to illustrate 
how such effects can influence the determination of the unitarity triangle. In Fig. 4, the 
constraint to the vertex of the unitarity triangle is given from \Vub\, Attib and e. Here the 
new physics effect can change the value of \Vtd\ and the shape of the bounds from e. It is 
different from the ref. [^ in which \Vtd\ was fixed and taken the value \Vtd\ = 0.0084. In 
our numerical calculations, we take \C,c\ = 9.8 and |^(| = 0.54, the mass of charged Higgs is 

11 



fixed at M^ = 200 GeV. The relative phase between them is taken to be 7r/3 and 27r/3 as 
two examples. The other input parameters are: B^ = 0.75±0.15, [Kifel/lKbl = 0.08±0.02 
and \Vcb\ = 0.04. For a comparison, we carry out a similar calculation in the SM with the 
same parameters, the result is plotted in Fig.5. It is found that the shape of the triangle 
can be largely changed when taking different values of the relative phase between C,c and 
^f The angle (3 of the triangle may be extremely small when cos 6' is close to 1. 

The second case, ex receives contributions from both the charged-scalar and the 
neutral-scalar exchanges. This case is more important than the one in which only the 
neutral-scalar exchange is dominant. This is because the relative phase between the the 
two contributions can largely affect the determination of \Vtd\- To illustrate such phase 
effect, we choose the ratio between the charged- and neutral-scalar contributions to Aitib 
to be 2:1, and take four typical values for the relative phase between them as 0, vr/S, 27r/3 
and TT. As it was pointed out by Soares and Wolfenstein [^] that if such phase emerges. 



then the unitarity angle extracted from B — > J/i/jKs will be the total phase 0m rather 
than /3. Here (pM is defined by 

with the index 'SM' and 'NEW denoting the contributions from the SM and new physics. 
Recently, the preliminary measurement of -B ^ J'ipKs was reported by CDF Collabora- 
tion [^ with a^Ks = 0-7Q-0.U- However, the resulting constraints on new physics are 



not very strong ||35| due to the large errors. A more precise measurement is expected at 
B factories. 

In Fig. 6, we plot the value of Vtd extracted from Atub in the p — 1] plane without 
considering the uncertainty of Bk{Bk = 0.75). The four curves are corresponding to 
the above four cases. The figure shows that the additional phase from y/^ • can strongly 
change the value of Vtd- Its modulus varies in the interval between 0.7 and 1.2 in this 
situation. 

As an example, the influence on the determination of the unitarity triangle in the case 
3 (see fig. 6) is plotted in Fig. 7. The bounds from e become lower than that from the SM 
due to the relative phase effect. As the three bounds from \Vub\, ArriB and e still have area 
in common, the triangle remains to be closed. However, as we have mentioned above, 
when the angle /3 is extracted from B — > J/ipKs, its value will be the total phase 0m 
which may be much larger. As a consequence, it could make the unitarity triangle to be 
'open. This possibility has been shown to happen in the case 3, where tan 0m is three 
times as large as tan/3. 
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FIG. 4. The constraints on the unitarity triangle in r] — p plane, the two different triangles 
correspond to 6 = it/ 3 (case a) and 9 = 27r/3 (case b). Where 9 is the relative phase between 
^c and ^f Other parameters are Bk = 0.75 ± 0.15, iKbl/l^bl = 0.08 ± 0.02 and \Vcb\ = 0.04. 
The allowed region between the two dashed curves are from \Vcb\, the allowed region between 
two dotted curves are from ArriB and the allowed region between two solid curves are from |ex| 
(the two curves below are for the case a for 9 = tt/3 and the two curves above are for the case 
b for ^ = 27r/3.) 
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FIG. 6. Constraints on Vtd from Atub- The relative phase between charged- and neu- 
tral-scalar exchange is taken to be 0(case 1, solid), 7r/3( case 2, dashed), 27r/3(case 3, dash-dotted) 
and 7r(case 4, dotted) 
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FIG. 7. The constraints on the unitarity triangle in case 3, where tan(/>M is three times as 
large as tan/3 

IV. CONCLUSIONS 

In conclusion, we have studied one of the simplest extentions of the standard model 
with an extra Higgs doublet, which we have simply called as an S2HDM, some constraints 
on the parameters in the S2HDM have been obtained from F^ — F° mixing processes. It 
has been shown that in general ^u^ ^c^ ^t and \C,c\ ^ |Ct|. A much larger D^ — D^ mixing 
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than the SM prediction is possible. Various sources of CP violation have been discussed. 
Their influences on the determination of the unitarity triangle have also been studied in 
detail. We found that angle 13 of unitarity triangle could be largely suppressed due to 
the new contribution from Higgs box diagrams. The phase from neutral Higgs exchange 
could strongly affect the extraction of {3 from B -^ J/ipKs- In some cases, such an effect 
could be so large that the unitarity triangle cannot remain to be closed. In particular, it 
may even result in the angle (3 determined from fitting the quantities \Vub\i ArriK and e 
being different from the one extracted from the decay process B -^ J/ipKs- If it is so, a 
clear signal of new physics is indicated. 
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Appendix 

In this appendix, we present some functions and quantities appearing in the text. 
A. The integral functions from the box diagrams with W boson and charged scalar 
exchanges 

B \x,x) = Vxx { — I — - 



-4 2 1 -a; A{l-xyx-x' 

+ [x '^ X ) - -— -\ 

^ ' ^{\-x){\-x'Y 

tdWW I \, ^ i X 

B X = - H -: + --, TTT In a; 

^ ' 4 41 -x 2 l-a;2 2 1 - a; P 



,2 



■ In y + (y ^ y') + 



{y-y'){l-yY ^ - ^> (I - y)(X - y') 

l + y , 2y , 



(l-yy (l-yf 
By {y,y , yw) = levr m^(-/6 + m^h) - 



'4"' ^ ^^ (l - y)il - y')il - yw) 

y{yw-y/4:) ^ yw , , ,. 
■ In — + [y ^y) 



(y - y')i^ - y)iy - yw) y 



Bi^^iy,yw) = lQn'mj,i-h + mi^h 



^ ^ ' {l-y){y-yw) 

{I - y)'^{l - yw) A{yw -y)'^{l-yw) y 

Bs [y,y) = —i jrr, T^iny-(y^y)- 



{y-y'){l-y)^ " ^" "^ {I - y){l - y') 






Bl'^{y,y\yw) = lQ7i^mlh 



{l-y){l-y'){l-yw) 
15 



{y - y'){^ - y){y - yw) 



^ yw , I ,^ 
In \- [y ^y) 
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\ny 
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where the functions Jj {i = 1,- ■ ■ ,6) are the euchdean integrals |l36[| : 
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B. The quantities used for calculating the CP-violating parameter e. 
Bi,{mi,m^) = ^M'A\\^^?U^B^''iy^^yi) + '^R<ii^';)B^'^{y^,yJ,yw)} 
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with 

The involved functions arise from the box graphs. From gluonic penguin graph with 
charged scalar exchange, we have 



and 
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